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ABSTRACT 

We present absorption line strength maps of 48 representative elliptical and lenticular galaxies 
obtained as part of a survey of nearby galaxies using our custom-built integral-field spectro- 
graph, SAURON, operating on the William Herschel Telescope. Using high-quality spectra, 
spatially binned to a constant signal-to-noise, we measure four key age, metallicity and abun- 
dance ratio sensitive indices from the Lick/IDS system over a two-dimensional field extending 
up to approximately one effective radius. A discussion of calibrations and offsets is given, 
as well as a description of error estimation and nebular emission correction. We modify the 
classical Fe5270 index to define a new index, Fe5270s, which maximizes the usable spatial 
coverage of SAURON. Maps of Hf3, Fe5015, Mg b and Fe5270s are presented for each galaxy. 
We use the maps to compute average line strengths integrated over circular apertures of one- 
eighth effective radius, and compare the resulting relations of index versus velocity dispersion 
with previous long-slit work. The metal line strength maps show generally negative gradients 
with increasing radius roughly consistent with the morphology of the light profiles. Remark- 
able deviations from this general trend exist, particularly the Mg b isoindex contours appear to 
be flatter than the isophotes of the surface brightness for about 40% of our galaxies without 
significant dust features. Generally these galaxies exhibit significant rotation. We infer from 
this that the fast-rotating component features a higher metallicity and/or an increased Mg/Fe 
ratio as compared to the galaxy as a whole. The H/3 maps are typically fiat or show a mild pos- 
itive outwards radial gradient, while a few galaxies show strong central peaks and/or elevated 
overall H/3-strength likely connected to recent star-formation activity. For the most prominent 
post-starburst galaxies even the metal line strength maps show a reversed gradient. 

Key words: galaxies: bulges - galaxies: elliptical and lenticular, Cb - galaxies: evolution - 
galaxies: formation - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

We are carrying out a survey of the dynamics and stellar popu- 
lations of 72 representative nearby early-type galaxies and spiral 
bulges based on measurements of the two-dimensional kinemat- 
ics and line strengths of stars and gas with SAURON, a custom- 
built panoramic integral-field spectrograph for the William Her- 
schel Telescope, La Palma I Bac on et alJl200ll hereafter Paper I). 
The goals and objecti ves of the SAURON survey are described in 
Ide Zeeuw et alJ ^2002, hereafter Paper II), which also presents the 
definition of the sample. The full maps of the stellar kinematics 
for the 48 elliptical (E) and lenticular (SO) galaxies are given in 
lEmsellem et alJ 120041 hereafter Paper III). The morphology and 
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kinem atics of the ionised gas emission are presented in lSarzi et alJ 
fcood hereafter Paper V). The stellar and gaseous kinematics of 
the spiral bulges are described in Falcon-Bar roso et alj |2006). Here 
we present maps of the absorption line strength measurements for 
the 48 E and SO galaxies in the survey. The analysis of the line 
strength maps, and a full analysis of the spiral galaxies in the sam- 
ple will be presented in later papers of this series. The data and 
maps presented here will be made available via the SAURON WEB 
page jhttp : / / www . strw . leidenuniv . nl/ sauron/| 

The measurement of absorption line strengths in combina- 
tion with stellar population models has been used for many years 
to probe the luminosity-weighted age, metallicity and abundance 
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ratios of ce rtain elements in integrated stellar population s (e.g. , 
FabeJl973tfBurstein et alJl984tlRosdi985lBrodie & Hanes|l98ti 
Bica et alll990[lGonzalezll993l-|Davies et alll993tlWorfhevll9'94l : 
Vazdekis 1999; Thomas et al. 2003). Integral- field spectroscopy al- 
lows to obtain spectroscopic information over a contiguous area on 
the sky and thus to identify two-dimensional structures (see e.g . , 
lEmsellemet"ai]ll99d IPeletier et al1ll999t: Idel Burso et"alll2001l) . 
For the first time we can apply this technique to absorption line 
strength observations in a large representative sample of early-type 
galaxies. 

One of the most widely-used methods for investigating the 
optical spectra of old (age ^ 1 Gyr) integrated stellar popula- 
tions in early-type galaxie s is to derive line strength indices in 
the Lick/IDS sys tem jBurstein et ail Il984t Iworthev et alJ Il994l : 
iTrager et al. 1998), and compare them with stellar population mod- 
els. We adopt this system in order to allow comparison with existing 
data. For future analysis we envisage the use of higher resolution 
line strength systems once they become available. 

In this paper we present the SAURON line strength measure- 
ments for our 48 representative E and SO galaxies, together with 
a brief analysis which focuses on the overall characteristics of the 
maps. A full stellar population analysis of the line strength maps 
with the help of stellar population models will be presented in a 
future paper in this series. 

This paper is organized as follows. In Section|2|we summarize 
our observational campaign, while the data reduction steps are out- 
lined in Section [5] along with a discussion on the flux-calibration 
of the data. In Section [4] we describe the definition and measure- 
ment of the line strength indices. The two-dimensional line strength 
maps are presented in Section|5| while in Section|6|we discuss line 
strength gradients and aperture corrections. In Section^we discuss 
average line strength measurements extracted from circular aper- 
tures of one-eighth effective radius, respectively, and compare them 
to the literature. The conclusions follow in Section|8| We comment 
on individual galaxies in Appendix lAl 



2 OBSERVATIONS 

The SAURON survey of 72 galaxies was carried out during eight ob- 
serving runs over four years (56 observing nights allocated). Details 
on the observing conditions, instrument set-up and exposure times 
are given in Paper III. The SAURON sample of 48 elliptical (E) and 
lenticular (SO) galaxies is representative of nearby bright early-type 
galaxies (cz ^ 3000km s ; Mb ^ — 18mag). As discussed in 
Paper II, it contains 24 galaxies in each of the E and SO subclasses, 
equally divided between 'field' and 'cluster' objects (the latter de- 
fined as belonging to the Virgo cluster, the Coma I cloud, and the 
Leo I group), uniformly covering the plane of ellipticity e versus ab- 
solute blue magnitude Mb- Tables Al and A2 in Paper II provide 
basic information on the objects. 

The SAURON survey is carried out with the low resolution 
mode of the instrument, giving a field of view of 33" x 41", con- 
tiguously sampled by 1431 0'.'94 x 0'.'94 square lenses, each of 
which produces a spectrum. Another 146 lenses sample a small 
region 1.'9 from the field center, which are used for simultaneous 
sky subtraction. The spectral resolution of the observed wavelength 
range 4800-5380 A is ~4.2A (FWHM, a inst ~ 108 km s" 1 ) and 
sampled at 1.1 A pixel" 1 . This wavelength range includes a set of 
potential emission lines (e.g. H/3, [OIII], [Nl]) and stellar absorption 
lines (e.g., H/3, Mg, Fe) which can be used to investigate the stellar 
populations. 



Each galaxy field was typically exposed for 4 x 1800 s, 
dithered by about 1". In 30% of the cases, we constructed mosaics 
of two or three pointings to cover the galaxy out to about one ef- 
fective radius R c , or, for the largest objects, out to ~0.5J? C - The 
footprints of these pointings are shown overlaid on a Digital Sky 
Survey image in Figure 1 of Paper III. 

In order to allow for inter-run calibration and provide tem- 
plates for redshift and velocity dispersion measurements we ob- 
served during each run a number of stars covering a broad range of 
spectral types. Specifically, we included stars from the Lick stellar 
library catalogue <Worthevetalll994 in order to calibrate our line 
strength measurements to the Lick/IDS system and its associated 
models (e.g., Worthev 1994, see also Section l4.2.3> . Spectrophoto- 
metric standard stars were also observed to calibrate the response 
function of the system (see Table0. 



3 DATA REDUCTION 

We reduced the SAURON observations with the dedicated 
XSAURON software developed at CRAL-Observatoire, and de- 
scribed in Paper I. The basic reduction steps include bias and dark 
subtraction, extraction of the spectra using a fitted mask model, 
wavelength calibration, low frequency flat-fielding, cosmic-ray re- 
moval, homogenization of the spectral resolution, sky subtraction, 
and flux calibration. The wavelength calibration is accurate to 0.1 A 
(6 km s - , rms). 

Multiple exposures of the same galaxy were merged and mo- 
saiced together. In this process the wavelength domain was trun- 
cated to a common range and the datacubes re-sampled to a spatial 
scale of 0'.'8 x 0>.'8 with orientation North up and East to the left. In 
order to allow for a meaningful analysis the final data-cubes were 
spatially binned to a minimum S/N of 60 per A using an adaptive 
scheme developed by Capp ellari & Copir] 12003). In this approach, 
the spectra are co-added by starting from the highest S/N lenslet, 
and accreting additional lenslets closest to the current bin centroid. 
A new bin is started each time the target S/N is reached. The re- 
sulting bin centroids are then used as starting points for a centroidal 
Voronoi tessellation, ensuring compact non-overlapping bins and a 
uniform S/N in faint regions. The binning used in this paper is iden- 
tical to the one used in Paper III and V. 

In the following we describe in more detail how the flux cali- 
bration was derived, since this is an important issue for the deriva- 
tion of line strength indices which was not covered in detail in pre- 
vious papers of this series. We also summarize the removal of emis- 
sion lines which is another important step before the measurement 
of absorption line strength. 



3.1 Flux calibration 

3.1.1 Overall flux calibration 

During each observing run a number of spectrophotometric stan- 
dards were observed, typically close to the center of the field of 
view (hereafter FoV) of SAURON (see Table0. As a first step we 
evaluated the inter-run consistency of the spectral response function 
of SAURON. For this purpose we investigated the relative changes 
in the spectral response function in repeat observations of six flux 
standard stars. In this analysis we are only interested in the relative 
changes so we eliminated the overall throughput difference by nor- 
malizing each set of observations of the same star in a 100 A region 
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Figure 1. lcr standard deviation of the spectral response function between 
our seven observing runs as a function of wavelength. The 22 stellar obser- 
vations used for this analysis were normalized in a 100 A region at 5060 A. 
The evaluation of the scatter was performed in ~ 10 A steps. 

in the middle of the wavelength range. In total, 22 repeat observa- 
tions 1 were used. 

The results are presented in Figure Q The relative run- to-run 
deviations in spectral response are smaller than 1% (la standard 
deviation) over the full wavelength range. An analysis of the varia- 
tions in absolute throughput for 13 standard stars observed in pho- 
tometric nights during the survey gives a standard deviation of 5%. 
No significant run-to-run variations could be established. Since the 
overall stability of the instrument is satisfactory we adopt a com- 
mon flux calibration curve for all runs which is described in the 
following. 

In order to establish an absolute flux calibration for our spec- 
tra a two step procedure was adopted: i) a relative spectral re- 
sponse function was established which maps throughput variations 
on ~ 100 A scales introduced by a wavelength selection filter in the 
instrument (see Pap er I); ii) the flux calibration tables compiled by 
iBohlin et all feOOlt "CALSPEC" directory) were used to establish 
an absolute calibration of the SAURON system. 

The first step in the calibration procedure is necessary since 
small-scale variations cannot easily be removed by the coarse wave- 
length steps of classical flux calibration tables. Furthermore, it is not 
uncommon to find large residuals at the position of absorption fea- 
tures (e.g., H/3). To overcome these problems we observed in runs 
4, 5, & 7 the white dwarf EG 131 (see Paper III). This star shows 
an exceptionally smooth continuum over the SAURON wavelength 
range, and in particular does not show any notable H/3 absorption 
(M. Bessel, private communication). This makes it an ideal cali- 
bration star to remove small-scale continuum variations. Unfortu- 
nately, there is no absolute flux calibration available for EG 131, 
but the energy distribution of a black-body with a temperature of 
~1 1800 K can be used to define the continuum (M. Bessel, private 
communication). 

Figure|2|shows the result from our analysis of the flux standard 
stars. We find several "wiggles" with a peak-to-peak amplitude of 
up to 3%. These wiggles, primarily derived from the observations of 
EG 131, can be confirmed in an independent way with the standard 
star G191B2B since Bohlin et al. (2001) provide a model flux cal- 
ibration curve at a similar resolution as the SAURON observations. 
This model flux curve is of very high quality and provides accurate 

1 Note, that some flux standard stars were observed multiple times within 
one run. 
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Figure 2. (a) The flux calibration curve for the SAURON survey normalized 
to arbitrary units, (b) Small-scale variations of the flux calibration curve 
after removal of a linear fit. 



Table 1. List of spectrophotometric standard stars 



Name 


Run observed 


Calibration 


(1) 


(2) 


(3) 


Feige 34 


5,6,7 


STIS+Oke 


Feige 66 


1,3 


FOS+Oke 


G191B2B 


1,4,5,6 


Model 


HZ 44 


1,3 


STIS 


BD+28 4211 


2,4 


STIS 


BD+33 2642 


7 


FOS+Oke 


HD 09352 


5 


STIS 


EG 131 


4,5,7 





Notes: Listed are the spectrophotometric standard stars which were 
used to establish the flux calibration curve for the SAURON system. 
Column ( 1 ) lists the name of the star while column (2) gives the run 
where the star was observed (see Paper III). Column (3) denotes 
the type of calibration file which was used for the absolute flux 
calibration (see Bohlin et al. 2001, for details). The star EG 131 is 
a special continuum correction star. See text for details. 



flux values even at small wavelength intervals. Both standard stars, 
EG 131 and G191B2B give consistent results. In order to establish 
the absolute flux calibration curve for the SAURON survey we used 
all flux standa rd stars listed in Ta ble [T] and the corresponding flux 
tables given in Bohli n et alJl200lh . 

Before applying the overall flux calibration curve to our data 
the observed spectrum is corrected for airmass according to the av- 
erage extinction curve for La Palma l Kins 1985). The spectra are 
not corrected for Galactic extinction since the correction would be 
negligible over the observed wavelength range. 

Although we derive from the photometric standard star obser- 
vations a very stable behavior of the instrument and a well cali- 
brated flux curve, we cannot directly transfer this conclusion to the 
full data set. For example, variations as a function of field of view 
and sky-subtraction errors for low signal-to-noise regions can have 
a significant effect on the flux calibration accuracy. 
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3.1.2 Continuum correction for line strength indices 

An analysis of our H/3 line strength measurements showed that we 
find non-physically strong absorption strengths in the outer bins of 
some galaxies often at the right side (as shown in Figure [5} of the 
FoV. A closer inspection revealed that the abnormal H/3 absorption 
strength is mainly caused by a wrong continuum shape across the 
H/3 region. In the outer bins typically many individual lenslets need 
to be summed up in order to reach our minimum target S/N of 60. 
Therefore, any systematic errors, such as imperfect tracing of the 
spectra on the CCD or sky subtraction, will be present in the final 
binned spectrum. Unfortunately we were unable to remove these er- 
rors in our data reduction and thus decided to correct the continuum 
in the binned spectrum. 

In order to obtain a good template spectrum for each galaxy 
we derive the optimal template for an average spectrum, within a 
10" radius, of each galaxy. This optimal template is the weighted 
combination of 19 spectra, spanning a range in a ge and metallicity , 
from the library of stellar population models of Vazdekis 1 1999), 
and additionally six spectra from the JonesJ 1 19971) stellar library 
from which the Vazdekis models are built to provide spectra with 
large Mg absorption strengths (for details see Paper III). The op- 
timal template is obtained from an emission line corrected (see 
Section l3~2t average spectrum, by simultaneously fitting the stel- 
lar kinematics together with the templates, and using a low-order 
polynomial. 

For each individual bin of the galaxy in question we fit this 
optimal template, with the kinematics fixed to the values of Pa- 
per III, together with an II th order multiplicative polynomial con- 
tinuum over the full wavelength range of the spectrum. Then in a 
second step we check if the continuum fit shows significant second 
and higher order structure over the wavelength range of the H/3 in- 
dex. Linear continuum variation does not affect the measured line 
strengths. Only for bins covering more than one lenslet and where 
we detect a significant 2 continuum variation we divide by the con- 
tinuum fit and thus apply the continuum correction before measur- 
ing the H/3 index. All bins which have been continuum corrected 
show a flag in the data-cubes of the public release and furthermore, 
we added (in quadrature) a constant systematic error of 0. 1 A to the 
H/3 index of corrected bins. The typical fraction of corrected bins 
in a data-cube is 4% while the extremes go from 1% to 33%. As an 
example Figure|3|shows the full H/3 map of NGC 3379 before and 
after the continuum correction. 

In Figure [4] we show a simulated long-slit observation of the 
H/3 absorption strength along the major axis of NGC 3379 derived 
from the SAURON data before and after the continuum correction. 
The top panel clearly shows the points at large negative radii where 
the continuum has been corrected. The literature comparison (bot- 
tom panel) shows a good agreement. 

The success of our continuum correction is critically depen- 
dent on the ability of the optimal template to fit the observed spec- 
trum. In general this may be a problem if the galaxy shows non- 
solar abundance ratios which are not mimicked by the template 
library. However, in the case of H/3 the index is only weakly de- 
pendent on abundance ratio iTrager et all200CtlThomas et al. 2003 ) 
and thus the method works well. Continuum variations can in prin- 
ciple also affect other indices such as Fe5015 and Mgb. Maps of 
Fe5015 are clearly affected, while Mgfe maps show evidence of 
continuum problems for only a handful of galaxies. For metal in- 

2 The rms deviations are greater than 0.006 after the removal of a straight 
line fit. 
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Figure 3. The H/3 map of NGC 3379 before (top plot) and after (bottom 
plot) continuum correction. The range in H/3 strength is indicated by the 
color bar and the numbers to the top and bottom of the color bar give the 
range in H/3 [A]. The scale is the same for both plots. The yellow dots 
in the bottom plot indicate all bins which were corrected for the continuum 
variation. North is up and east to the left. The spatial axis are given in arcsec. 

dices abundance ratios are a critical issue. The continuum vari- 
ations and unaccounted abundance ratio differences between the 
galaxy spectrum and the optimal template are largely degenerate 
(at least for our data and models). Thus we decided not to correct 
the Fe5015, Mg b and Fe5270s indices. For selected galaxies we re- 
move bins which suffer from significant continuum problems (see 
also Section|5J. 

3.2 Emission-line corrections 

Many early-type galaxies show signs of relatively weak emis- 
sion from i onized gas and dust in the interstellar medium (e.g., 
Goudf rooii et alJll994h . Three out of four of the main absorption 
line features in the SAURON wavelength range are potentially af- 
fected by emission lines; namely, H/3, [Oni]A4959, [OIII]A5007, 
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Figure 4. Effects of continuum correction for NGC3379 (major axis). The 
top plot shows SAURON data before (open circles) and after (filled circles) 
continuum correction for the H/3 index. The bottom panel shows a compar- 
ison with lite rature data. Error b ars are only plotted for the SAURON data. 
The data of iDavies et all h993j) is shown as crosses and mirrored about 
the zero point; data from lGonzalezl (19931 and private comm unication) is 
show n as open triangles ; the d ata of McDermid et al. 1 2005) is shown as 
stars. |Sanchez-Blazrjuez 1 2004 ) observations are shown as open diamonds. 
Literature data outside ±27" is not shown. 

and [NI]A5200. Before measuring absorption line indices these 
emission lines, where present, need to be removed taking into ac- 
count that the gas kinematics can be decoupled from the stellar kine- 
matics. This is done as part of our efforts to analyze the emission 
lines themselves. The emission-line measurements and results were 
described in detail in Paper V. Here we give only a short summary. 

For five galaxies we find no emission above our detection limit 
(NGC 821, NGC 2695, NGC 4387, NGC 4564, and NGC 5308); for 
a further seven galaxies we find only weak evidence for emis- 
sion (NGC 4270, NGC 4382, NGC 4458, NGC 4473, NGC 4621, 
NGC 4660, and NGC 5845), while all of the remaining galaxies 
in our survey show clear signs of emission. All emission lines are 
modeled with a simple Gaussian profile. Using the stellar kinemat- 
ics derived in Paper III we re-derive the best fitting, optimal stellar 
template 3 in combination with the best fitting emission template 
featuring [OIII] only. While the [OIII] strengths and kinematics are 

3 The composition of the template library wa s improved w ith respect to 
the one used in Paper III. Three stars from the IjoneJ Q 997) library were 
exchanged with more suitable ones which allowed a better match of the 
spectral characteristics in large early-type galaxies and thus an improved 
emission correction. The effects on the derivation of the stellar kinematics 
are negligible. 



free parameters in this first step, the regions of potential H/3 and 
[Nl] emission are masked. If no significant [OIII] emission is found 
we assume that there is also no H/3 and [Nl]. This limits our ability 
to measure weak emission associated with star-formation regions 
where H/3 / [OIII] ^> 1. This is an unavoidable consequence of 
the fact that weak H/3 emission cannot always reliably be measured 
without constraining its kinematics. Whether emission is found or 
not is determined by a measure of the amplitude of the emission 
lines with respect to the level of both statistical and systematic de- 
viations from the fit to the stellar continuum. As a cutoff we use an 
amplitude-over-noise (hereafter A/N) ratio of 4. This corresponds 
to a detection limit of 0. 1 A in the central, high S/N, parts of the 
line strengths maps and is reduced to approximately 0.2 A in the 
outer parts. 

If we find significant [OIII] emission we re-fit the galaxy spec- 
trum allowing for potential H/3 and [Nl] emission. For galaxies 
with weak emission lines we constrain the H/3 and [Nl] emis- 
sion kinematics to be the same as the [OIII] ones. For ten galax- 
ies (NGC 2768, NGC 3032, NGC 3414, NGC 4278, NGC 4374, 
NGC 4459, NGC 4486, NGC 4526, NGC 5838, and NGC 5846) 
with relatively strong emission we are able to fit the H/3 emission 
kinematics independently while the [Nl] kinematics are tied to H/3 
(see Paper V for details). The A/N detection limits for H/3 is set to 
3 and 5 for the constrained and unconstrained kinematics, respec- 
tively. The detection limits for H/3 emission are 0.06 A and 0.2 A in 
the inner and outer parts of the galaxies, respectively. 

For [Nl] the A/N detection limit is always set to 4. We note 
that especially for the [Nl] emission line the dominating source of 
noise is not shot noise from the spectrum, but template mismatch 
due to unaccounted for non-solar abundance ratios in the galaxies. 
We therefore measure the A/N ratio for the [Nl] emission line in the 
Mg b region which is a good guide for the degree of template mis- 
match. Due to the difficulty of quantifying the effects of template 
mismatch, the detection limits of the [Nl] line are not well estab- 
lished. They will be, however, similar to the [OIII] detection limits 
and somewhat larger for objects with large non-solar abundance ra- 
tios. 

In Figure|5|we show an example map for the emission correc- 
tion of the H/3 feature. NGC 4526 has an inclined dust disk in the 
central 20" x 20" (see Paper V) showing very regular gas kinemat- 
ics and a circumnuclear region with particularly strong H/3 emis- 
sion. After emission correction a relatively thin elongated region of 
strong H/3 absorption coinciding with the dust disk appears. This 
example is representative of the strongest emission corrections in 
our data. 

Many previous authors corrected for H/3 emission using an 
average ratio between the [OIII]A5007 and H/3 lines. For instance 
Trage r et all <200fJ) concluded that 0.6 times the [OIII]A5007 emis- 
sion is a reasonable estimate for the H/3 emission. However, in their 
sample of 27 galaxies the correction factor varies from 0.33 to 1.25 
and it is therefore doubtful whether this correctio n is accurate for 
an individual galaxy (see also M ehlert et all200(3) . We confirm this 
with our data (see Figure |6j and find a median correction factor 
of 0.67 with a range from 0.3 to 2.3. There is weak evidence that 
the average correction factor is a function of emission strength, be- 
coming smaller at large [OIII]A5007 strength. More importantly, 
we point out that emission line ratios generally vary significantly 
across the galaxies (see error bar in Figure [6] and [OlII]/H/3 maps 
in Paper V). NGC 3032 shows an exceptionally large ratio of H/3/ 
[OIII] ~ 2.1, pointing towards ongoing star-formation over a sig- 
nificant area of the galaxy. 



© 0000 RAS, MNRAS 000, 000-000 



6 Kuntschner et al. 



NGC4528 




NGC4536 H)S emission corrected [OIII] A5007 [A] 




■30 o 20 



Figure 5. The H/3 map of NGC 4526 before (top plot) and after (bottom 
plot) emission and continuum correction. The range in H/3 strength is indi- 
cated by the color bar and the numbers to the top and bottom of the color bar 
give the range in H/3 [A]. The scale is the same for both plots. North is up 
and east to the left. The spatial axis are given in arcsec. The small circular 
region of strong H/3 absorption visible in both maps (Ra ~ 18" and Dec ~ 
7") corresponds to the position of a star in the field of view. 

4 THE LINE STRENGTH MEASUREMENTS 

Here we describe which line strengths indices we use for the 
SAURON survey, the calibration of our measurements to the stan- 
dard Lick/IDS system and the evaluation of the line strength errors. 



4.1 Line strength indices 

In the Lick/IDS system absorption line strengths are measured by 
indices, where a central feature bandpass is flanked to the blue and 
red by pseudo-continuum band-passes (see Figure 0. The mean 
height in each of the two pseudo-continuum regions is determined 
on either side of the feature bandpass, and a straight line is drawn 
through the midpoint of each one. The difference in flux between 



Figure 6. Median H/3 emission vs [OIII] emission for all galaxies with sig- 
nificant detection of emission lines. The error bar in the top left corner shows 
a typical range of emission strength found within individual galaxies (esti- 
mated as a robust 1<t standard deviation). The dashed line reflects the rela- 
tion of H/3 = 0.67 X [OIII]A5007 corresponding to our median correction 
factor. The position of NGC 3032 is indicated by its name; for details see 
text. 

this line and the observed spectrum within the feature bandpass de- 
termines the index iTrager et al. 1998). For most absorption fea- 
tures, the indices are expressed in angstroms of equivalent widths 4 . 

The exact wavelength definitions of the indices, the observed 
wavelength range and the redshift of the object determine which 
indices can be measured. The full wavelength range of SAURON 
is 4760-5400 A. However, the tilting of the wavelength selection 
filter in order to avoid ghost images, reduces the usable wavelength 
range common to all lenslets to approximately 4825-5275 A (see 
Paper I). For the nearby sample of SAURON galaxies (i.e., cz ^ 
3000 km s _1 ), this wavelength range allows us to measure three 
Lick/IDS indices: H/3, Fe5015 and Mg b (see Tabled 

The usable wavelength range is a well-determined function of 
position on the CCD. The red wavelength limit is most affected, and 
varies from 5275 to 5380 A over the FoV. For lenslets with a red- 
extended wavelength coverage we can measure a further Lick/IDS 
index: Fe5270. It is desirable to measure this index since it is an 
important ingredient to investigate the Mg- to-Fe abu ndance ratios 
in integrated stellar populations (Wort hev et alll992T) . 

Since the red wavelength cutoff is a monotonic function across 
the field of view, by making the red bandpass of the Fe5270 index as 
blue as possible, we effectively increase the usable field of view. For 
this reason we have designed a new index for the Fe5270 absorption 
feature with a reduced red wavelength coverage. The bandpass def- 
initions of all indices used for the SAURON survey are summarized 
in Tabled 

The main aim of the re-definition of the Fe5270 index is to 
increase the effective field of view. A useful re-definition of the in- 

4 For broad molecular bands, the index is expressed in magnitudes. No such 
indices are used in the SAURON system. 
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Figure 7. The set of line strength indices shown on a central SAURON spec- 
trum of NGC 3379. The shaded regions are the central band-passes whereas 
the dashed lines indicate the side-band-passes. The solid lines show the 
pseudo-continuum determined from the side-band-passes. The location of 
the [OIII]A5007 emission line is indicated by the dotted line. 

dex should also, however, obey the following conditions: (i) small 
and well-determined line of-sight-velocity-distribution (hereafter 
LOSVD) corrections (see Section l4,2,2t : (ii) good Poisson statis- 
tics; (iii) similar sensitivity towards changes in age, metallicity and 
abundance ratios as the original Lick/IDS Fe5270 index in order to 
allow easy conversion. 

In order to evaluate the optimal new index definition we first 
created a catalogue of model spectra resembling typical SAURON 
observations of early-type galaxies. In fact we employed the same 
model spectra library as used for the optimal template determina- 
tion (see Paper III) and the LOSVD corrections (see Section l4.2.2t . 
By using this model library we cover reasonably well the range of 
expected observations for the Fe5270 absorption feature. 

All six wavelength values (see Table |2j which define an in- 
dex were systematically changed in 0.5 A steps until an optimal 
definition" 1 was found. The bluest wavelength definition of the red 
pseudo-continuum which still gives reasonable LOSVD corrections 
and error statistics for the index turns out to be 5308.0 A. This is a 
reduction of « 10 A in the red-most limit of the original Lick/IDS 
definition, equivalent to an increase in the field of view of 13%. 

The LOSVD correction for a purely Gaussian LOSVD for the 
new index is only 5.1% at 200 km s _1 as compared to 10.3% for the 
original index. While the index could be improved in terms of its ve- 
locity dispersion corrections and Poisson statistics it also shows a 
good linear relation with the original Lick/IDS index, i.e., it shows 
a similar dependence on age and metallicity (see Figure Fortu- 
nately, the re-definition of this index did not significantly change its 
overall sensitivity to abundance ratios (Kom et al. 2005). The best 
fitting linear relation determined from stellar population models of 
Vazdekis 1 1999) with ages > 4Gyr is: 

Fe5270 = 1.26 x Fg5270 s + 0.06. (1) 

The formal errors of the fit are smaller than ±0.01. In this paper 
we use the above equation to convert our Fe5270s measurements to 

5 The index definition was optimized at the instrumental resolution of 
SAURON (kj 108 km s~ 1 ) in order to allow for the new generation of stellar 
population models which will be able to make predictions at this spectral 
resolution. 




Fe52?0s (SAURON) [A] 

Figure 8. The relation between the new Fe5270g index and the original 
Lick/IDS index Fe5270. The small plus signs represent SAURON obser- 
vations of NGC 3379 within an eight arcsec radius of the center, while the 
open diamonds show the central measure ments of a c omplete sample of 
early-type galaxies in the Fornax cluster I Kuntschner 2000). The thick solid 
red line is a linear fit to the SAURON da ta: Fe5270 = 1.2 8 X Fe5270 s + 
0.03. Overplotted are model predictions by Vazdekis 1 1999) at the Lick/IDS 
resolution. 

the Lick Fe5270 index. We note, an empirically determined conver- 
sion formula compares well (Fe5270 = 1.28 x Fe5270s + 0.03; see 
Figure[8}. 

In order to maximize the information for the line strength in- 
dices we produce for each galaxy two data-cubes: one containing 
the H/3, Fe5015 and Mg b indices for the full field of view and one 
presenting the new Fe5270s index for a reduced field of view. Due 
to the adaptive two-dimensional binning, the final spatial sampling 
of the two data-cubes is not identical. 



4.2 Calibration to the Lick/IDS system 

In order to allow a meaningful comparison between stellar popula- 
tion model predictions and observed data, the index measurements 
need to be carefully calibrated to the Lick/IDS system. There are 
three effects to account for: (a) the difference in the spectral reso- 
lution between the Lick/IDS system and the SAURON instrumental 
set-up; (b) the internal velocity broadening of the observed objects 
and (c) small systematic offsets caused by continuum shape differ- 
ences iWorthev & Ottavianll99l : lKuntschnerl2000l) . 

In the following paragraphs we describe the individual steps 
taken to ensure an accurate calibration of the SAURON data to the 
Lick/IDS system. 

4.2.1 Spectral resolution correction 

First we need to adjust for the difference in spectral resolution be- 
tween the Lick/IDS system and the SAURON observations for line 
strength measurements. The nominal resolution of the Lick/IDS 
system in the SAURON wavelength range is 8 — 9 A (FWHM) with 
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Table 2. Bandpass definitions of SAURON line strength indices 



Index Blue pseudo-continuum Central bandpass Red pseudo-continuum Units Source 



H/3 


4827.875 


-4847.875 


4847.875 


- 4876.625 


4876.625 


-4891.625 


A 


Tfasrer et al. ( 19981 


Fe5015 


4946.500 


-4977.750 


4977.750 


- 5054.000 


5054.000 


-5065.250 


A 


Trageretal. ( 19981 


Mgb 


5142.625 


-5161.375 


5160.125 


-5192.625 


5191.375 


- 5206.375 


A 


Trageretal. (19981 


Fe5270 


5233.150 


-5248.150 


5245.650 


- 5285.650 


5285.650 


-5318.150 


A 


Trageretal. ( 19981 


Fe5270 s 


5233.000 


- 5250.000 


5256.500 


- 5278.500 


5285.500 


- 5308.000 


A 


This paper 



a mean of 8.4 A iWorthev & Ottavianill997() . We degrade our spec- 
tra to approximately the Lick/IDS resolution by broadening to an 
instrumental resolution of a — 212 km s _1 . 

4.2.2 Line of -sight-velocity -distribution corrections 

The observed spectrum of a galaxy is the convolution of the in- 
tegrated spectrum of its stellar population(s) by the instrumental 
broadening and line of-sight velocity distribution (LOSVD) of the 
stars. These effects broaden the spectral features, in general reduc- 
ing the observed line strength compared to the intrinsic values. In 
order to compare the index measurements for galaxies with model 
predictions we calibrate the indices to zero velocity dispersion and 
the nominal Lick/IDS resolution. 

In the past, most authors have only taken into account the first 
moments, v and a, of the LOSVD to correct the line strengths in- 
dices. However, with the availability of high quality measurements 
of the higher order terms hi and hi which describe non-Gaussian 
deviations (see e.g., Paper III) we need to apply more accurate cor- 
rections. Typically a non-zero /13 term does not affect the Lick in- 
dices significantl y, but a non-zero hi term can have significant ef- 
fects lHaJl998LKuntschnej2504Tl . 

In order to determine the corrections for each index we use the 
optimal template determined for each spectrum individually during 
the derivation of the kinematics (for details see Paper III). Gener- 
ally, this optimal template is a good representation of the galaxy 
spectrum. For each spectrum a correction factor, Cj(a, /13, hi), is 
determined such that 

CjipMM) =I j ((T = 0,h 3 =0,hA = 0)/I j (a,h 3 ,h4,), (2) 

where Ij is the index measured from the optimal template con- 
volved to the Lick resolution and additionally convolved with the 
LOSVD given in brackets. A LOSVD corrected index is then 

jcorr = C] (a,h 3 ,h 4 ) X Ij aW . 

4.2.3 Lick/IDS offsets 

Although we matched the spectral resolution of the Lick sys- 
tem, small systematic offsets of the indices introduced by contin- 
uum shape differences are generally present (note that the origi- 
nal Lick/IDS spectra are not flux calibrated). To establish these off- 
sets we compared our measurements for star s in common with the 
Lick/IDS stellar library iWorthev et al. 1994). In total we observed 
73 different Lick/IDS stars with 174 repeat observations 6 . 

Figure [5] shows the difference between Lick/IDS measure- 
ments and ours, where repeat observations of the same star were 
averaged. The mean offsets and associated errors for each index are 



evaluated by a biweight estimator to minimize the influence of out- 
liers and are summarized in Table|5| The formal error in the offset 
is taken to be the standard error on the mean a I \/N s tars- 

As an independent test of this procedure, we compare the 
Lick/IDS measurements for the 42 galaxies in common with our 
survey 7 . This is more difficult than the similar comparison for stars 
since for galaxies, e.g., aperture differences and seeing variations 
can affect the results. In order to match the Lick/IDS standard aper- 
ture of l'.'4x4" we extracted the luminosity-weighted mean of a cir- 
cular aperture with radius r — l'.'335 from each line strength map. 
It would in principle be better to extract the exact Lick aperture: this 
is however not possible since the position angles of the Lick obser- 
vations are unknown to us. For this comparison the SAURON mea- 
surements were not corrected for emission but internal line broad- 
ening was taken into account (see Section l4.2.2t . 

The results are overplotted in Figure [5] as open diamonds and 
the offsets, evaluated with a biweight estimator, are listed in Table|5| 
The overall offsets inferred from the galaxies are in good agreement 
with the offsets derived from the stars. To further test the offsets 
for Lick indices derived from flux calibrated data we made use of 
the ljonesl ^997) library. First the data was broadened to the spec- 
tral resolution of the Lick/IDS system with a wavelength dependent 
Gaussian assuming a constant spectral resolution of FWHM = 1.8 A 
for the Jones stars. Then we compared the index measurements for 
the 128 stars in comm on between the Jones library and the Lick 
observations iWorthev et ai]ll994ll . The offsets and associated er- 
rors, derived with a biweight estimator, are listed in Table [3] col- 
umn 4 (see also. lNorris et alj2006l) and are in good agreement with 
the offsets derived from the SAURON observations for the H/3 and 
Fe5015 indices. The agreement for the Mgb and Fe5270 indices 
is less good, but differences are still small with offsets < 0.1 A. 
Furthermore, our determination of Lick offsets derived from Jones 
stars i s in excellen t agreement with an earlier investigation carried 
out bv lWorthev & Ottavianl 1 19971 Table 9). Although the offsets 
are small, they can significantly change the age and metallicity esti- 
mates. Particularly, the H/3 index presents a problem, since one can 
find with the currently available stellar population models age esti- 
mates older than the age of the universe if offsets are applied. In the 
end we decided to apply Lick offsets only for the H/3 and Fe5015 
indices (see Table|3| column 5) since all three offset determination 
methods give consistent results. For the Mg b and Fe5270 indices no 
offsets are applied since the SAURON stellar observations indicate 
offsets consistent with zero. For all indices we quote conservative 
offset errors which reflect the information derived from the stars 
and galaxies. 



Since the stars have only small relative velocities and were typically ob- 
served at the center of the SAURON field of view we were able to measure 
the original Lick/IDS index Fe5270 for all stars. 



7 In the Lick/IDS survey the Mg b, Fe5015, and H/3 indices are not avail- 
able for one, three and two galaxies, respectively. 
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Figure 9. Comparison of Lick/IDS and our measurements for 73 stars in common (plus signs). For one star (HD 148783, spectral type M6III) our measurements 
were drastically different for all four indices. The dashed line shows the mean offset derived by a biweight estimator (see also Tablel3*l. Observations of galaxies 
in common between Lick and SAURON are shown as open diamonds. A typical error bar for an individual star is shown in the upper left corner of each panel. 
For the galaxy comparisons we plot a separate representative error bar. The error is dominated by the observational error of the Lick/IDS data. Note that the 
stellar observations span very well the index range covered by the galaxies. 



Table 3. Offset estimates to the Lick/IDS system 



Index 


Offset (stars) 


Offset (galaxies) 


Offset (Jones stars) 


Applied offsets 


(1) 


(2) 


(3) 


(4) 


(5) 


H/3 


-0.13 ±0.03 A 


-0.06 ± 0.04A 


-0.12 ±0.02 A 


-0.13 ±0.05 A 


Fe5015 


+0.28 ±0.05 A 


±0.27 ±0.09 A 


±0.23 ±0.04 A 


±0.28 ±0.05 A 


Mgb 


-0.01 ±0.03 A 


±0.01 ± 0.04A 


-0.08 ±0.02 A 


±0.00 ±0.05 A 


Fe5270 


±0.02 ±0.02 A 




-0.07 ±0.02 A 


±0.00 ±0.05 A 



Notes: Column (1) gives the index name, while column (2) and (3) give the mean offset (Lick - 
SAURON) to the Lick/IDS system evaluated from the stars and galax ies in commo n, respectively. 
Column (4) shows the offsets derived from the comparison of the Ijonesl 119971) library to the 
Lick/IDS observations for stars in common. Column (5) shows the final offsets we applied to our 
data and their mean adopted error. 



4.3 Estimating errors for line strength indices 

Several authors have e valuated the sources of errors for line strength 
indices ICardiel et alll 19981 and references therein). They can be 
divided into random and systematic errors. The treatment of sys- 
tematic errors is generally very difficult, while the random errors 
can be determined quite accurately. In the following paragraphs we 
concentrate on the random errors. 

The main ingredient to evaluate the random errors for line 



strength indices in the SAURON system are the noise spectra which 
are provided by the data reduction procedure for each individual 
lenslet (see Paper I). The noise spectrum represents the errors prop- 
agated through the full data-reduction until the final, merged and 
binned data-cube from which velocities, velocity dispersions, line 
strength indices etc. can be determined. 

Before measuring a given index on a spectrum one needs to 
know: a) the average recession velocity in order to place the index 
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band-passes correctly and b) a measure of the LOSVD in order to 
correct the index for broadening effects (see Section 14.2.21 . Both 
parameters have associated uncertainties (for details see Paper III) 
which will contribute to the final index error. 

For each galaxy, the random error for each index measurement 
was estimated via a Monte Carlo approach. For every spatial bin, N 
realizations of the associated spectrum were produced based on its 
corresponding noise distribution (assumed to be Gaussian). Like- 
wise, N values of the recession velocity and velocity dispersion 
were chosen within the uncertainties of their measured values, thus 
creating a total number of N 3 realizations of the measured index, 
from which the rms uncertainty was determined. A typical number 
for N in our analysis was 30. For bins where emission in H/3, [OIII] 
or [Nl] was detected and subtracted we add the associated errors in 
quadrature (see Section l3~2l and Paper V). 

Concerning the systematic errors we consider two main con- 
tributions: (i) errors from the continuum correction for the H/3 line 
strengths and (ii) the error in the calibration to the Lick system. The 
continuum correction errors (see Section 13.1.21 assumed to be at 
a constant level of 0. 1 A) are added to the affected bins while the 
Lick system offsets are quoted in Table [5] and are not added to the 
individual data points, but shown as global offset errors in figures. 



5 OBSERVED LINE STRENGTH MAPS 

Figures [T7all71l below present maps of the absorption line strengths 
of the 48 objects, ordered by increasing NGC number. For each 
galaxy we show the total intensity reconstructed from the full 
wavelength range of the SAURON spectra (see also Paper III), and 
the two-dimensional line strength distributions of Mg b, Fe5015, 
Fe5270s and H/3, overplotted with isophotes of the reconstructed 
image spaced by single magnitude steps. The maps are all plot- 
ted with the same spatial scale, and oriented with respect to the 
SAURON field for presentation purposes. The relative directions of 
north and east are indicated by the orientation arrow next to the 
galaxy title (we note that the orientation and binning of the maps is 
identical to Papers III and V). The maximum and minimum of the 
plotting range is given in the tab attached to each parameter map, 
and the color bar indicates the color table used. In order to allow 
for an easy comparison between galaxies the plotting range of the 
H/3 maps is fixed to 1.0 - 2.9 A (with the exception of NGC 3032, 
3156, 4150 which show larger H/3 absorption strengths in the cen- 
tral parts). The colors are adjusted such that blue to green shades 
correspond to the H/3 strength predicted for old stellar popula- 
tions (~12Gyr), while stronger H/3, corresponding to the presence 
of younger stellar populations, are r epresented by red and yellow 
shades (see e.g.. lThomas et alJl2003h . For the metal line maps we 
use an independent plotting range for each galaxy in order to better 
visualize the line strengths gradients across the maps. 

For some galaxies we found non-physically large values of the 
Fe5015 line strength in the outer regions where many individual 
lenslets are averaged to achieve the target S/N of 60. The same ef- 
fect can be seen to a much lesser extent in Mg b maps. Similar to the 
continuum variations which cause some of the H/3 measurements to 
be corrupted in the outskirts, we suspect that the Mg b and Fe5015 
indices are also affected (see Section |3. 1 .21 . Due to the degener- 
acy between abundance ratio variations and continuum effects we 
cannot apply a continuum correction for the metal lines. Therefore, 
we decided to remove the most affected bins by hand from the fi- 
nal data-cubes. These bins are indicated with grey color in the line 
strength maps presented in Figures fl7all71l 



5.1 Overview of the line strength maps 

The maps in Figures [1 7al 1 711 show a wealth of structures and we 
give comments on individual galaxies in Appendix^ Some general 
trends are apparent and we discuss these briefly in the following. A 
full stellar population analysis of the line strength maps with the 
help of stellar population models will be presented in a future paper 
in this series. 

The metal line strength maps often show negative gradients 
with increasing radius roughly consistent with the morphology 
of the light profiles. Remarkable deviations from this trend exist: 
NGC 3032, 3156, 4150, and 4382 show a central depression in Mg b 
line strength. This structure is always accompanied by a strong peak 
in the H/3 maps indicating the presence of a recent post-starburst or 
even low level ongoing star-formation. Interestingly, the Fe maps 
of post-starburst galaxies do not always show a depression such as 
the Mg 6 index which can be understood in the reduced sensitivity 
to stellar population age of Fe indices. 

Significant, but low-level, deviations between isoindex con- 
tours of the metal indices and the isophotes are discussed sep- 
arately in Section 15.21 We find that enhanced Mg b absorption 
strength is connected to fast-rotating components (as defined in 
Capp ellari et all2006l hereafter Paper IV) in the galaxies. Perhaps 
the strongest case in our sample is presented in the Mgb map 
of NGC 4570 where the Mg b strength along the major axis runs 
through a dip outside the bulge region and then begins to rise to- 
wards larger radii where the galaxy kinematics start to be dominated 
by an outer disk (see also iFisher et all 1 996. for a long-slit study of 
SOs). 

For galaxies with weak H/3 absorption (~ 1.4 A) the maps 
show typically mild positive gradients or are consistent with being 
flat. A few galaxies (e.g., NGC 3489, NGC 7332 and NGC 7457) 
show H/3 maps with an overall elevated H/3 line strength indicat- 
ing a spatially extended, recent star-formation episode. Other galax- 
ies (NGC 4382, NGC 4459 and NGC 4526) exhibit a more centrally 
concentrated region of relatively strong H/3 absorption pointing to- 
wards a more localized star-formation event. 

Regions of strong H/3 absorption are often associated with dust 
features and detected in emission (see emission maps and unsharp- 
masked images in Paper V). The case of NGC 4526 deserves a spe- 
cial note, since here one can see the direct connection between an 
almost perfect dusty disc seen close to edge on (apparent in the 
reconstructed image), and the location of strong H/3 absorption in- 
dicating the presence of young (~1 Gyr) stars. 



5.2 Mg b isoindex contours versus isophotes 

One of the most interesting aspects of integral-field spectroscopy is 
the capability to identify two-dimensional structures. For the first 
time we can use this in connection with line strength indices and 
compare isoindex contours with the isophotal shape. One might ex- 
pect the index to follow the light in slowly-rotating giant elliptical 
galaxies, since the stars are dynamically well mixed. However, one 
could imagine that dynamical substructures with significantly dif- 
ferent stellar populations could leave a signature in line strength 
maps which are sensitive to e.g., metallicity. 

The Mg6 index is the best determined index in our survey 
and potential differences between isoindex contours and isophotes 
should be most apparent. Indeed we find a number of galaxies 
where the isoindex contours clearly do not follow the isophotes, 
e.g., NGC 4570 and NGC 4660. Both galaxies show fast-rotating 
components (see Paper III) along the direction of enhanced Mg b 
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Figure 10. Top row: The interpolated Mg b map and the reconstructed im- 
age of NGC 3377. Bottom row: The best fitting elliptical model with con- 
stant position angle and ellipticity to the Mg b map and the reconstructed 
image. The best fitting ellipticity for the isophotes and the Mgb map is 
0.473 ± 0.003 and 0.573 ± 0.022, respectively. The x- and y-axis are given 
in arcsec; North is up and East to the left. 



strength. In order to further study such a possible connection we de- 
termine the best fitting (x 2 sense) simple, elliptical model for each 
of the reconstructed images and Mg b maps in our survey. For this 
model we impose a constant ellipticity and position angle as func- 
tion of radius, where the position angle is determined as the average 
position angle of the reconstructed image. An example of this pro- 
cedure is shown in Figure fTol We then compare the derived charac- 
teristic ellipticities for the isophotes and Mg b maps with each other. 
Clearly such a simple model does not give a good fit to all galaxies, 
but it provides us with a robust, global shape measurement of the 
isoindex contours as compared to the isophotes. 

In order to avoid confusion by the presence of significant 
dust absorption we exclude the following galaxies from the anal- 
ysis: NGC 2685, NGC 3032, NGC 3156, NGC 3489, NGC 4374, 
and NGC 4526. Furthermore we exclude NGC 4486 because of the 
presence of the non-thermal emission from the well-known jet. This 
leaves a sample of 41 galaxies. 

The results of the ellipse fitting are shown in Figure fTTI Errors 
are evaluated via 50 Monte Carlo simulations per map where the 
Mg b strength of each bin is varied according to the index errors. 
Additionally, we allowed for a 5% fraction of the data values to be 
considered as outliers and varied those data points by 5 times its 
original error. We required a minimum error of 0.08 A per bin 8 and 
also excluded all bins covering more than 20 individual lenses from 
the analysis. 



8 In the central, high S/N regions of the Mg b maps our internal error es- 
timates are very small and we use a minimum error of 0.08 A to reflect 
systematic errors. 




0.2 0.4 0.6 0.8 

Mean ellipticity (photometry) 

Figure 11. Comparison of the average ellipticity of constant Mg b strength 
with the best fitting elliptical model of the isophotes. Errors are evaluated by 
a Monte-Carlo simulation and represent lcr errors. All galaxies which are 
more than 2a above the one-to-one line are indicated by filled circles and 
their galaxy names. 



Applying a 2a error cut we find 16 out of the 41 galaxies 
to be consistent with no deviation between isophotes and isoin- 
dex contours. Seven galaxies appear to have rounder isoindex con- 
tours than the isophotal shape (NGC 2695, 2699, 2768, 3384, 3414, 
4262, and 5845). However, 18 galaxies appear to have more flat- 
tened Mgb contours than the isophotes (NGC 821, 2974, 3377, 
3608, 4278, 4382, 4473, 4477, 4546, 4550, 4564, 4570, 4621, 4660, 
5831, 5838, 5982, and 7332). Most of these galaxies show a high 
degree of rotational support and significant /13 terms. Thus, the 
flattened Mg b distribution suggests that the fast-rotating compo- 
nents in these galaxies exhibit a stellar population different from 
the main body. The enhanced Mg b strength can be interpreted to 
first order as higher metallicity and/or increased [Mg/Fe] ratio. Fur- 
thermore, it is evident from the line strength maps that the detailed 
isoindex shapes cannot always be modeled with a simple elliptical 
structure. For example, in NGC 4570 one can see a central bulge 
component with a strong radial Mg b gradient while at larger radii 
the onset of a disk component w ith strong Mgb is visible (e.g., 
Ivan den Bosch & Emse llem 1998). A detailed analysis of the line 
strength distribution and its possible connection to the kinematics 
will be presented in a forthcoming paper. 



6 LINE STRENGTH GRADIENTS AND APERTURE 
CORRECTIONS 

6.1 Line strength gradients 

A number of previous studies have investigated line str ength 
gradi e nts with the help of long- slit data (e.g.. iGorgas et alj 
| l99ri ICarollo et alJll993t iDavies et aljll993l : iFisher et alJll996t 
iMehlert et alj|2003l) . Line strength gradients can be used to study 
the formation history of early-type galaxies since different forma- 
tion models predict different gradients. In a nutshell, monolithic 
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collapse models <Carlberd fl984) predict steep metallicity gradi- 
ents with metal rich centers whereas hierarchical models, follow- 
ing a merger tree, predict shallower gradients due to the dilution of 
any line strength gradients existing in the pre-merger units (White 
1980). 

We now derive line strength gradients from the SAURON sam- 
ple. First, all indices are expressed in magnitudes (as is the well 
known Mg2 index) which is indicated by a prime sign [']. The con- 
version between an index measured in A and magnitudes is 



NGC3379 



index' = —2.5 log ( 1 — 



index 



AA J 



) 



(3) 



where AA is the width of the index bandpass. 

In Section l5~2l we have mentioned the fact that overall, the 
morphology of the line strength maps resembles the corresponding 
photometry. This motivated our choice to derive robust and simple 
line strength gradients by averaging the indices along lines of con- 
stant surface brightness (isophotes) with equal steps in log flux. The 
average line strength in each radial bin is derived after applying a 
3<r clipping algorithm. The radius for each bin is calculated as the 
median major axis radius normalized to the effective radius along 
the major axis a c — Re/yl — e, where R c is the effective radius 
and e is the average ellipticity of the galaxy within the SAURON 
field. The effective radii R c (see Table[5J are derived with a R 1 ^ 
growth curve analysis from our wide-field MDM (1.3m) imaging 
survey of the E/SO sample, supplemented by archival HST images. 
The photometry and analysis thereof will be presented in a forth- 
coming paper of this series. 

For each galaxy we fit an error-weighted straight line to the 
index values at a given radius such that gradients are defined as 



A index' = 



5 index ' 
S log(R/R c ) ' 



(4) 



The error of the index value in each surface brightness bin is taken 
as the lcr scatter of the data-points within this bin after applying a 
3<7 clipping algorithm. We further restrict the fitting range to radii 
between 2" and R c . An example fit to the data of NGC3379 is 
shown in Figure IT2l 

In Figure [T3l we present for the full sample the fitted gradi- 
ents for the H/3, Fe5015, and Mg b indices plotted against the av- 
erage (luminosity-weighted) index values within a circular aperture 
of one eights of an effective radius (Re/8, see Section Q for de- 
tails). For elliptical galaxies the H/3 gradients have a mean value of 
0.000 ± 0.006, i.e., they are consistent with a flat relation (see also 
Meh lertet al] |2003). The lenticular galaxies show a range in H/3 
gradients where stronger (negative) gradients are correlated with in- 
creasing central (R c /8) H/3 absorption strengths. This is mostly a 
consequence of the presence of galaxies in our sample which harbor 
young stars in the central regions, but have little or no star forma- 
tion in the outer parts. If these young stars amount only to a small 
fraction in mass, the gradients will within a few Gyrs become shal- 
lower and the central H/3 absorption strength will decrease, moving 
the galaxies back to the bulk of points in Figure fOl 

We find mostly negative gradients for the metal line gradients, 
consistent with the literature (e.g.. | Davies et alll993tlCarollo et al] 
ll993tlKobavashi & Arimotdll9991) . Median gradients for Fe5015' 
and Mgb' are -0.014 and -0.025, respectively. Using the full sam- 
ple we find significant correlations between the line strength gra- 
dients and the average values within (R c /8) for Fe5015 and Mgfe. 
The correlation coefficient derived from a (non-parametric) Spear- 
man rank-order test is -0.55(< 1%) and -0.49(< 1%), for the 
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Figure 12. Line strength gradients for NGC 3379 as function of log (R/R c ) 
along the major axis. Error bars reflect the lcr scatter of data-points along a 
given band of constant surface brightness. The red line represents the best- 
fitting straight line to the data. The radial extent of the red line reflects the 
fitting range in log R/R c . The fitted line strength gradients are noted in the 
lower left of each panel. 



Fe5015' and Mg&' indices, respectively. The probability that the 
parameters are not correlated is given in parentheses. 

A possible correlation between the Mg gradient (as measured 
by the Mg2 inde x) and the central M g line strength had been sug- 
gested bv lGonzalez & Gorgasl il99 5'. see also Carollo et al. 1993), 
but some other investigations failed to detect a significant correla- 
tion (e.g.. iKobavashi & Arimotdll999l : iMehlert et alj f2003). New 
long-slit observations of a sample of 82 early-type galaxies pre- 
sented by Sanchez-Blazquez (2004) do confirm the trend found by 
iGonzalez & GorgasUl995ll 

In our sample the correlation for the Mg6' index is largely 
driven by the (lenticular) galaxies with Mg b' < 0.13. These galax- 
ies are the ones with young stars in the central regions. This can be 
seen in Figure fT4b where we show the relation between H/3' and 
Mg b' gradients. There is also a good correlation between the gradi- 
ents of_fhe_MgJ/andFe50 15' indices as shown in Figure fl4k (see 
also lFisheretalJ l996). This indicates that both indices are sensitive 
to similar stellar population parameters (e.g., metallicity). 

In summary it appears that in our sample there are signif- 
icant correlations of the line strength gradients with the central 
line strength in early-type galaxies. The relations are predominantly 
driven by the lenticular galaxies and furthermore by the presence of 
young stellar population in the central regions of galaxies. Environ- 
mental differences do not seem to play a major role. It is now inter- 
esting to explore if the above correlations also hold if line strength 
gradients are compared to index strength at larger radii, i.e. more 
representative regions for the galaxies as a whole. 

In Figure fT3l (right hand side panels) we plot the fitted gradi- 
ents for the H/3, Fe5015, and Mgb indices against the index value 
at 1/4 effective radius along the major axis. The index value at 1/4 
effective radius is taken to be the value predicted by the fit. We have 
chosen this particular reference point, since it is within the range of 
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Figure 13. Left column of panels: Index gradients as function of luminosity-weighted index values within a circular aperture of one eights of an effective 
radius (R c /8). The red filled and open circles are cluster and field ellipticals, respectively; blue filled and open triangles are cluster and field SOs, respectively. 
The smaller symbols indicate galaxies for which our data cover less than half the effective radius. Error bars are formal errors from the straight line fits. A 
gradient of 0.0, i.e. a fiat relation, is indicated by the dashed lines. Median gradients are given in the lower left of each panel. See text for details. Right column 
of panels: Index gradients as function of (fitted) line strength at 1/4 effective radius along the major axis. Symbols are the same in both diagrams. 



the observed line strength gradients for all galaxies and is more rep- 
resentative of the galaxy as a whole rather than being biased to the 
central regions. 

In this new diagram we do not find a significant correlation 
for the Fe5015' index. For the H/3 index there is still a correlation 
visible but it is exclusively driven by the three galaxies in the lower 
right part of the plot. The situation for Mg b' appears more complex. 
There is no significant correlation for the full sample while the SOs 



show a trend similar to the one derived for the 1/8 effective radius 
values. The elliptical galaxies even appear to follow a positive cor- 
relation between line strength gradient and Mg b' strength at R e /4. 
However, the latter trend is not significant. 

We conclude that for early-type galaxies in our sample there is 
a correlation between line strength gradients of H/3, Fe5015 & Mg b 
and centrally averaged line strength values (R c /8). These correla- 
tions are driven to a large extent by the presence of young stellar 
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Figure 14. (a) Fe5015' gradients versus Mg b' gradients, (b) H/3' gradients 
versus Mg b' gradients. The red filled and open circles are cluster and field 
ellipticals, respectively; blue filled and open triangles are cluster and field 
SOs, respectively. The smaller symbols indicate galaxies for which we do not 
cover half the effective radius. Error bars are formal errors from the straight 
line fits. 



populations in the central regions of (typically lenticular) early-type 
galaxies and thus are caused by galaxies which have recently ex- 
perienced secondary star-formation which may account only for a 
minor fraction of the total galaxy mass. 



6.2 Aperture corrections 

For studies of galaxies over a significant redshift and size range 
aperture corrections have to be applied in order to allow for a fair 
comparison of kinem atic and stellar population parameters (e.g., 
J0reen ~n et alj|l995l) . In Paper IV we present an aperture correc- 
tion for the luminosity weighted second moment of the line of-sight 
velocity distribution (a) while here we present aperture corrections 
for the line strength indices H/3, Fe5015, and Mg b. 

For this we show in Figure [131 the luminosity-weighted line 
strength measurements within an aperture of radius R, normalized 
to their value at half of the effective radius. The profiles were mea- 
sured by co-adding the SAURON spectra within circular apertures of 
increasing radius for all 39 early-type galaxies which have data out 
to at least half of the effective radius. If the SAURON FoV does not 
cover the complete circular aperture at a given radius, we calculate 
an equivalent radius of a circular aperture which matches the area 
of the SAURON FoV. 

The metal line indices show generally a decreasing profile with 
increasing (R/R c ), while the H/3 index shows a wide range in 
aperture corrections. Typical aperture corrections are small where 
a correction from 1/4 effective radius to one effective radius would 
amount to approximately 0.5% for the H/3 index and ~9% for the 
metal indices. However, it is clear that substantial galaxy-to-galaxy 
variations in the profiles exist. Particularly, for galaxies harboring a 
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Figure 15. Luminosity-weighted index measurements within an aperture of 
radius R, normalized to its value at 1/2 R c . Each galaxy with data out to 
at least 1/2 R c is shown with the thin colored lines. Data points inside 
R = 2" are excluded to minimize seeing effects. The black line is the me- 
dian power-law relation (indexn/index noTIn ) = (R/R n oim) p where p 
is given in Table l4l The robustly estimated la errors of the power law fits 
are indicated by the dashed black lines. See text for details. The three galax- 
ies with (Mg 6 R /Mg 6 e/2 ) < 1.0 at small radii are NGC 3032, NGC 3156 
and NGC 4150. 



young stellar population in the center (i.e. positive Mg b gradients) 
the application of an average aperture correction can be severely 
wrong. 

The best fitting power-law relations to the line strength mea- 



surements as function of normalized radius 



were deter- 



mined for each galaxy individually by fitting all available data out- 
side R — 2". Then a median power-law index and a robustly es- 
timated scatter was derived. These median power-law indices for 
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Table 4. Aperture corrections for line strength indices 



Table 6. Linear fits to index-cr relations (for log cr 2.1) 



norm ) = (R/R n 
(Fe5015 R /Fe5015 norm 

) = (R/R 

n 

(Mgft R /Mg 

"norm 

) = (R/R 



\+0.001±0.029 
i) 

\-0.067±0.025 



-0.061±0.028 



each absorption feature are only provided as a convenient approx- 
imation of the true aperture corrections and presented in Table |4| 
They are in reasonable agreement with previously published aper- 
ture corrections tforgensen et all 1 995: Meh lert et all2003h . 



7 INTEGRATED INDEX-er RELATIONS 

After the presentation of the complete line strength maps (Sec- 
tion |5j and the average gradients (Section |6| we now present 
luminosity-weighted line strength indices averaged over a circular 
aperture. We focus on one key diagnostic: the index-cr relations. 

Our integrated measurements were derived by averaging the 
(luminosity-weighted) spectra within a circular aperture with a ra- 
dius of (iZe/8), but imposing a minimum aperture of 2'.'4 x 2'.'4 
(3x3 pixel). From this high S/N spectrum we re-measure the 
kinem atics (V, cr, and h± ) as described in Paper III (see also 
ICappellari & Emselleml2004h and then determine the line strength 
as presented in Section|4] Because of limited field coverage we do 
not quote central averaged Fe5270s measurements for eleven galax- 
ies (NGC474, 524, 821, 2695, 2768, 3608, 4270, 5198, 5308, 5846, 
and 5982). In a further eight galaxies the (R c /8) aperture is not 
fully covered by our maps, however, only a few bins, correspond- 
ing to less than 10% coverage, are missing and we judge the average 
spectrum to be still a good representation of the galaxy. 

In order to provide a more global measurement for each galaxy 
we also derived a central averaged spectrum from all data available 
within one effective radius. Because of the sometimes severely lim- 
ited field coverage of the Fe5270s index we do not quote one effec- 
tive radius averages for it. Furthermore, since our line strength maps 
of the other indices do not always cover the full area of one effec- 
tive radius we need to asses the degree of necessary aperture cor- 
rections. The galaxies with the smallest coverage (NGC 4486 and 
NGC5846) feature line strength data out to about 30% of R c . Even 
in these extreme cases we estimate corrections to be small (< 8%; 
see Section l6~2l . For the 32 galaxies for which the data extends to 
less than one effective radius, we apply the corrections of the form 
given in Table|4] The median coverage of the line strength maps is 
0.8-Rc- A subset of one R c index measurements is used in Paper IV 
to determine stellar M / L ratios and compare them with high fi- 
delity dynamical estimates of the global M/L inside one R a . 

The final, luminosity-weighted line strength measurements 
within a circular aperture of (R c /8), and within R c , are given in 
Table|5| The formal line strength errors are smaller than ±0.01 A 
for all galaxies. We note, however, that there are systematic errors 
which we estimate to be of the order of 0.06, 0.15, 0.08, 0.06 A for 
the H/3, Fe5015, Mg6, and Fe5270s indices, respectively. For the 
luminosity weighted averages of the velocity dispersion we adopt 
an error of 5% (see Paper IV). 

Our main result for the (R c /8) extractions is shown in Fig- 
ure ^| For the index-cr relations all indices are expressed in 
magnitudes (see Section loTY The Mgb index shows a tight cor- 



H/3' 


= -0.103(±0.017) logo- 


+ 0.301(±0.039) 


Mg b' 


= +0.153(±0.021)logcr 


- 0.195(±0.050) 


Fe5015' 


= +0.000(±0.010) logo- 


+ 0.080(±0.024) 


Fe5270' 


= +0.002(±0.010) logtr 


+ 0.075(±0.024) 



Notes: The linear fitting was performed by taking into account 
errors in both variables. The error in the intercepts and slopes 
have been derived by scaling the errors in the observables by a 
constant factor until the \ 2 probability is 0.5. 



(e.g., Terlevich et alJl98ltlBurstein et alll988: Bender et alll993l: 


Colless et alJ 


1999 




Bernardi et al. 


2003 




Worthev & Collobert 



120031: iDenicolo etafl l2005h . The H/3 index shows a relatively 
tight, negative correlation for galaxies with velocity dispersions 
log o ^ 2.1. For lower velocity dispersions there is evidence for 
an inc reased scatter (see also e.g.. lKuntschneJ200filCaldwell et all 
2003). The deviations from the main relation for galaxies with 
log cr < 2.1 clearly anti-correlate in the Mg b, H/3-cr relations. 

The two Fe indices show roughly constant index values for 
velocity dispersions log a ^ 2.1. Similar to the Mg b and H/3-cr re- 
lations we also find an anti-correlation with respect to the fitted rela- 
tion between Fe and H/3 index strength at a given cr albeit at lower 
significance. At lower velocity dispersions we find there is weak 
evidence for an increasing index strength with increasing velocity 
disper sion. Similar trends have been noticed bv lProctor & Sansoml 
12002) in an analysis of spiral bulges and E/SOs. The absence of 
a clear correlation between Fe indices and the ce ntral velocity dis- 
persion has been noted by several authors (e.g.. lFisher et aljll996l : 
J0rge nsenlll999l) and we can confirm this for our sample. How- 
ever, for other samples a significant correlation between Fe indices 
and the central velocity dispersion has been foun d (e.g.. lKuntschnerl 
l200CtlCaldwell et alj2003tlBemardi et all2003l) . 

In order to compare our data with the literature we choose 
two datasets which cover luminous early-type galaxies lGonzale3 
Il993t iTrager et all |2000L (Re/8) extractions) and fainter ones 
1 Caldwell et al. 2003, global extractions). Our sample of 48 early- 
type galaxies shows good agreement with the literature for the H/3, 
Mg b and Fe5270 ind ices (see F igure 1 161 . For Fe5015, which was 
only observed by Caldwell et al. (2003), an offset of approximately 
0.01 mag at intermediate velocity dispersions is visible, in the sense 
that our data has the larger values. The best fitting linear index-cr 
relations for log cr ^ 2.1, using errors in both variables are given in 
Tabled 

There is no significant evidence for an environmental influ- 
ence on the index-cr relations in our sample. However, it is gener- 
ally the low velocity dispersion SOs which show the strongest anti- 
correlation between H/3 and metal line strengths. The reason for the 
scatter in the index-cr relations is typically ascribed to a combina- 
tion of age, metallicity and abundance ratio varia tions at a given a 
(e.g.,[Kuntschner et al. 2001; Thomas et al. 2005). We explore this 
further in a forthcoming paper in this series. 



8 CONCLUDING REMARKS 

The maps presented in this paper are the result of a comprehen- 
sive survey of the absorption line strength distributions of nearby 
early-type galaxies with an integral-field spectrograph. This data set 
demonstrates that many nearby early-type galaxies display a signifi- 
cant and varied structure in their line strength properties. This struc- 
ture is sometimes very apparent as in the case of the post-starburst 
affected central regions of galaxies, or it can be of subtle nature as 
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Table 5. List of line strength measurements within a circular aperture of R c /8 and R c 



Name 


R R 


ti a, x I -* **6 


c/8 


H/3 


Fe5015 


Mgb 


Fe5270 s 


H/3 


Fe5015 


Mg6 




1 J 




[km s _1 ] 


[Al 


[Al 


[Al 


rAi 


[Al 


[Al 


[Al 








-Rc/8 


Rc/8 


-Rc/8 




-Rc/8 


Re 


Re 


Rc 


f 1 ) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(11) 


NGC474 


29 


0.71 


163 


1.68 


5.45 


4.11 


- 


1.81 


4.85 


3.45 


NGC524 


51 


0.61 


249 


1.39 


5.57 


4.55 


- 


1.50 


5.38 


4.20 


NGC821 


39 


0.62 


202 


1.53 


5.32 


4.35 


- 


1.57 


4.65 


3.72 


NGC1023 


48 


0.56 


206 


1.49 


5.60 


4.51 


2.34 


1.51 


4.98 


4.19 


NGC2549 


20 


0.90 


142 


2.07 


5.85 


4.05 


2.40 


2.02 


5.10 


3.54 


NGC2685 


20 


1.34 


83 


2.01 


5.05 


3.57 


2.19 


2.04 


4.21 


3.04 


NGC2695 


21 


0.95 


223 


1.31 


5.05 


4.57 


- 


1.27 


4.34 


3.94 


NGC2699 


14 


1.41 


146 


1.83 


5.86 


4.28 


2.40 


1.76 


4.76 


3.58 


NGC2768 


71 


0.39 


205 


1.70 


5.02 


4.08 


- 


1.64 


4.43 


3.64 


NGC2974 


24 


1.04 


242 


1.64 


5.27 


4.40 


2.25 


1.71 


5.14 


4.09 


NGC3032 


17 


1.19 


100 


4.57 


4.58 


1.73 


1.62 


3.98 


4.02 


2.14 


NGC3156 


25 


0.80 


61 


3.92 


3.80 


1.60 


1.60 


3.00 


3.66 


1.79 


NGC3377 


38 


0.52 


142 


1.85 


4.95 


3.83 


2.09 


1.86 


4.37 


3.20 


NGC3379 


42 


0.64 


215 


1.43 


5.23 


4.57 


2.20 


1.44 


4.81 


4.15 


NGC3384 


27 


0.75 


158 


1.94 


5.92 


4.16 


2.47 


1.84 


5.12 


3.70 


NGC3414 


33 


0.60 


236 


1.44 


5.08 


4.48 


2.15 


1.50 


4.49 


3.72 


NGC3489 


19 


1.05 


100 


2.82 


4.99 


2.84 


2.01 


2.60 


4.40 


2.60 


NGC3608 


41 


0.49 


189 


1.49 


5.30 


4.36 


- 


1.60 


4.73 


3.75 


NGC4150 


15 


1.39 


78 


3.47 


4.23 


2.17 


1.69 


2.87 


4.11 


2.35 


NGC4262 


10 


2.06 


199 


1.42 


5.07 


4.54 


2.09 


1.49 


4.46 


3.90 


NGC4270 


18 


1.09 


137 


1.77 


5.09 


3.45 


- 


1.74 


4.67 


3.24 


NGC4278 


32 


0.74 


248 


1.23 


4.78 


4.67 


2.04 


1.41 


4.59 


4.21 


NGC4374 


71 


0.43 


294 


1.40 


5.19 


4.52 


2.22 


1.42 


4.66 


4.08 


NGC4382 


67 


0.38 


179 


2.16 


5.14 


3.45 


2.12 


1.98 


4.72 


3.36 


NGC4387 


17 


1.16 


96 


1.60 


5.03 


3.92 


2.18 


1.51 


4.42 


3.67 


NGC4458 


27 


0.74 


95 


1.63 


4.49 


3.73 


1.90 


1.59 


3.97 


3.27 


NGC4459 


38 


0.71 


180 


2.16 


5.39 


3.91 


2.27 


1.85 


4.68 


3.43 


NGC4473 


27 


0.92 


198 


1.49 


5.55 


4.59 


2.33 


1.50 


4.88 


4.03 


NGC4477 


47 


0.43 


168 


1.65 


5.15 


4.17 


2.21 


1.59 


4.71 


3.84 


NGC4486 


105 


0.29 


317 


1.14 


5.30 


5.12 


2.22 


1.15 


4.70 


4.62 


NGC4526 


40 


0.66 


240 


1.86 


5.59 


4.32 


2.29 


1.62 


4.93 


4.16 


NGC4546 


22 


0.94 


223 


1.54 


5.49 


4.59 


2.27 


1.54 


4.64 


3.92 


NGC4550 


14 


1.28 


81 


2.14 


4.57 


3.15 


2.00 


1.99 


4.32 


2.93 


NGC4552 


32 


0.63 


277 


1.39 


5.93 


5.04 


2.31 


1.37 


5.33 


4.55 


NGC4564 


21 


1.02 


173 


1.52 


5.98 


4.79 


2.42 


1.55 


4.85 


4.00 


NGC4570 


14 


1.43 


199 


1.45 


5.87 


4.70 


2.34 


1.45 


4.79 


3.99 


NGC4621 


46 


0.56 


229 


1.40 


5.53 


4.73 


2.28 


1.43 


4.75 


4.14 


NGC4660 


11 


1.83 


229 


1.43 


5.87 


4.88 


2.34 


1.47 


4.80 


4.07 


NGC5198 


25 


0.80 


209 


1.49 


5.50 


4.66 


- 


1.56 


4.70 


3.93 


NGC5308 


10 


2.04 


252 


1.46 


5.42 


4.57 




1.46 


4.92 


4.23 


JNOCjol 5 


52 


0.39 


15 1 


1.51 


j.lb 


4.35 


Tic 

2. ID 


1.52 


4.74 


1 H 

4. 12 


NGC5831 


35 


0.67 


166 


1.77 


5.55 


4.09 


2.35 


1.78 


4.78 


3.32 


NGC5838 


23 


0.87 


290 


1.64 


5.86 


4.55 


2.35 


1.58 


5.05 


4.19 


NGC5845 


5 


4.45 


285 


1.51 


6.12 


4.80 


2.38 


1.52 


5.38 


4.41 


NGC5846 


81 


0.29 


237 


1.34 


5.50 


4.85 




1.32 


4.89 


4.38 


NGC5982 


27 


0.94 


269 


1.63 


5.96 


4.50 




1.60 


5.18 


3.98 


NGC7332 


11 


1.91 


135 


2.28 


5.93 


3.71 


2.30 


2.16 


4.80 


3.17 


NGC7457 


65 


0.39 


62 


2.28 


4.81 


2.88 


2.02 


2.25 


4.33 


2.68 



Notes: (1) NGC number. (2) Effective (half-light) radius R c measured from HST/WFPC2 + MDM images. (3) Ratio 
between the maximum radius Rmax sampled by the SAURON observations and R c . The maximum radius is calculated 
as the equivalent circular radius matching the area of the full SAURON field of view. Note, that for five galaxies the values 
given here are slightly smaller than the ones given in Paper IV since we have removed bad bins for the line-strength 
indices (see also Section |5}. (4) Velocity dispersion of the luminosity weighted spectrum within a circular aperture of 
Re/8. (5) - (8) Fully corrected line strength index measurements of the luminosity weighted spectrum within .Rc/8 for 
the H/3, Fe5015, Mgb and Fe5270g indices. Due to limited field coverage we cannot determine the Fe5270g indices 
for eleven galaxies. (9) - (11) Fully corrected line strength index measurements of the luminosity weighted spectrum 
within -R c for the H/3, Fe5015, and Mgfe indices. For galaxies with less than one R c coverage we applied the aperture 
corrections given in Tablel4l Formal errors of the line strength indices are below 0. 1 A. We note, however, that there are 
systematic errors which we estimate to be of the order of 0.06, 0. 15, 0.08, 0.06 A for the H/3, Fe5015, Mg b, and Fe5270s 
indices, respectively. For the velocity dispersion we adopt an error of 5%. 
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Figure 16. (i? c /8) aperture line strength index versus a relations. All in- 
dices are given in magnitudes. The Fe5270s index was converted to the 
Fe5270 index by using Equation Q The filled and open circles are cluster 
and field ellipticals, respectively; filled and open triangles are cluster and 
field SOs, respectively. Mean error bars reflecting mostly systematic errors 
are shown in the upper left corner of each panel. Small dots are taken from 
Caldwell et al. 12003). The small plus signs represent the data of Gonzalez 
1 1993). Data from the literature with a < 50km s _1 are not shown. The 
dashed lines are fits to the data for all galaxies with log a ^ 2.1. The results 
of the linear fits are summarized in TablelSI 



seen in the deviations of Mg b isoindex contours compared to the 
isophotes of galaxies. 

The two-dimensional coverage of the line strengths allows us 
to connect the stellar populations with the kinematical structure of 
the galaxies and thus improve our knowledge of the star-formation 
and assembly history of early-type galaxies. For example, we find 
the Mgb isoindex contours to be flatter than the isophotes for 18 
galaxies which generally exhibit significant rotation. We infer from 
this that the rotational supported sub-structure features a higher 
metallicity and/or an increased Mg/Fe ratio as compared to the 
galaxy as a whole. Further steps in this direction will be presented 
in forthcoming papers of this series. 

The metal line strength maps show generally negative gradi- 
ents with increasing radius, while the H/3 maps are typically flat or 



show a mild positive outwards radial gradient. A few galaxies show 
strong central peaks and/or elevated overall H/3-strength likely con- 
nected to recent star-formation activity. For the most prominent 
post-starburst galaxies even the metal line strength maps show a re- 
versed gradient. We use the maps to compute average line strengths 
integrated over circular apertures of one-eighth effective radius. The 
resulting index versus velocity dispersion relations compare well 
with previous long-slit work. 
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Figure 17a. Maps of the line strength the 48 E and SO galaxies in the SAURON representative sampl e. The SAURON spectra have been spatially binned to a 
minimum signal-to-noise of 60 by means of the centroidal Voronoi tessellation algorithm o f lCappellari & Copinll20o3) . All maps have the same spatial scale. 
From top to bottom: i) reconstructed total intensity; ii) Mg b line strength; iii) Fe5015 line strength iv) Fe5270s line strength and v) H/3 line strength. 
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Figure 17b. 
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Figure 17c. 
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Figure 17d. 
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Figure 17e. 
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Figure 17g. 
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Figure 17h. 
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Figure 17i. 
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Figure 17j. 
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Figure 17k. 
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Figure 171. 
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APPENDIX A: DESCRIPTION FOR INDIVIDUAL 
GALAXIES 

Here we briefly comment on the line strength structures ob- 
served in the SAURON maps of the E/S0 sample presented in 
this paper. Individual comments on the stellar kinematics and 
the emission line maps are presented in Papers III and V, re- 
spectively. A thorough and quantitative assessment of the line 
strength maps with the help of stellar population models will 
be carried out in subsequent papers. Independently derived line 
strength maps for NGC524, NGC 1023 NGC3379, NGC3384, 
NGC4550, NGC 7332 and NGC 7457 based on our obse r vation s 
were presented a nd dis c ussed in[Afanasi ev & Sirchenkc] |2002|). 
Sil'ch enko et alJ (2002), Si l'chenko et alJ 1(2003) and Sil'chenko 

NGC 474: This galaxy (Arp 227), famous for its shell structures 
(e.g. lTurnbull et alll999h . shows normal metal line strength gradi- 
ents. The H/3 map appears relatively flat in the inner regions while 
we observe rising H/3 strength towards larger radii. 

NGC 524: This galaxy shows regular Mg b gradients and the H/3 
map shows a mild positive gradient. The two Fe maps show a noisy 
structure which may partly be caused by low S/N in the outer bins. 

NGC 821: A close to edge-on galaxy with a rapidly rotating disk 
like component (see Paper III) shows an enhanced core region in 
all metal line maps. There is also evidence for the Mg b contours to 
be flatter than the isophotes. The H/3 map appears rather featureless 
(see also lMcDermid et all2005t) . Th e maps are consistent with the 
results from long-slit spectroscopy of|Sanchez- BlazaueJ iiool . 

NGC 1023: This SB0 galaxy, with a prominent twist in the central 
velocity field (see Paper III), shows a central concentration in all 
metal lines. The metal line strength at larger radii appear flat and 
even tend to rise again at the edges of the field of view. The H/3 
map is consistent with a constant value over the full field of view. 

NGC 2549: This galaxy with a thin, rapidly rotating component 
(see Paper III) shows enhanced Fe5015 and Fe5270s strength along 
this component. The H/3 line strength is also elevated in the central 
regions. 

NGC 268 5: This f a mous object, the Helix galaxy 
iBurbidse & Burbidgd Il959t IPeletier & Christodouloij 1 19931) . 
shows significant extinction due to polar dust lanes on its North 
East side. The H/3 map shows enhanced values over large regions 
of the FoV, while the metal line strength maps are elevated along 
the rotation direction. 



NGC 2695: The metal line strength maps show normal gradients 
consistent with the isophotes. The H/3 map is flat, however, at no- 
tably low values. 

NGC 2699: This galaxy shows regular metal line strength gradi- 
ents consistent with the isophotes. The H/3 map is relatively flat at 
a slightly elevated level compared to the bulk of the objects in our 
sample. 

NGC 2768: This galaxy with a rather cylindrical velocity field 
shows regular metal line strength gradients. The H/3 map appears 
rather flat. Significant dust extinction is present north of the center 
<MicharJl999T) . 

NGC 2974: This rapidly rotating galaxy shows evidence of en- 
hanced metal line strength along the rotation direction. The 
H/3 map appea r s rela tively flat. The long-slit observations of 
iRampazzo et all l2005t) sh ow declining H/3 strength with radius 
while lCarollo et alj /l993) data suggest a modest increase with ra- 
dius. Our emission line maps (Paper V) show significant [OIII] and 
H/3 emission basically over the full FoV and thus different H/3 emis- 
sion corrections are likely the cause of the literature disagreement 
on H/3 absorption strengths. 

NGC 3032: This dusty galaxy shows negative Mg b and Fe5270s 
line strength gradients. The Fe5015 map shows a ring-like structure 
of weak absorption while the center exhibits strong Fe5015 absorp- 
tion. The H/3 map reveals a prominent peak in the center and shows 
evidence of rising values at larger radii. The unusual structure of 
the line strength maps and the relatively strong emission detected 
in this galaxy (see Figure[6|and Paper V) suggest that a recent star- 
burst has occurred in this galaxy with weak evidence of ongoing 
star formation. 

NGC 3156: Another dusty galaxy with negative Mgb line 
strength gradients in the center. The Fe5015 and Fe5270s maps do 
not show a clear structure. The H/3 map exhibits a strong peak in 
the center with a gradual decline to the outskirts. The line strength 
maps suggest a pos t-starburst scenario for this galaxy (see also 
iDenicolo et all2005h . 

NGC 3377: This dis ky galaxy jBender et alJ Il988l: 
Micha rd"& Simierl ^88) shows metal line strengths which 
are more flattened than the isophotes. Conversely, the H/3 absorp- 
tion has a rather constant moderate va lue across the field. T he line 
strength maps are consistent with Sanchez-Blazquez (2004). 

NGC 3379 (M105): This rather round object shows normal metal 
line strengths gradients consistent with the isophotes. The H/3 map 
shows weakly increasing line strength towards larger radii, as found 
by other authors (e.g., Davies et al. 1993; Sanchez-Blazquez 2004, 
see also Figure|4j. 

NGC 3384: This SB0 galaxy shows a dynamically cold compo- 
nent in the central 5" (see Paper III), which may correspond to th e 
red major axis disk component found by iBusarello et alJ ll996l) . 
The line strength indices were already presented in Paper II. The 
metal line strengths show regular gradients consistent with the 
isophotes. The H/3 absorption map, at moderately elevated values 
(see also lKuntschner et alj|200lh . appears flat in the inner regions 
with weak evidence for increasing H/3 strength at large radii. 

NGC 3414: This disturbed galaxy (Arp 162) contains a kinemat- 
ically decoupled component (see Paper III) within the central 10". 
The line strength maps appear to have regular gradients. 

NGC 3489: This dusty galaxy shows a fast-rotating central com- 
ponent and a complex morphology of the odd higher-order kine- 
matic moment /13 (see Paper III). The H/3 ab s orption st rength 
is glo bally high (see also lKuntschner et alJl200ltlRampazzo et alJ 
2005), with a sharp central (2 ") enhancement, also evident with 
OASIS jMcDermid et alj2004l) . 
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NGC3608: This galaxy shows a kinematically decoupled core 
(Paper III and references therein) extending to around 13". Our 
metal line strengths appear more flattened than the isophotes within 
this region. The H/3 absorption map shows a mild positive gradi- 
ent wi th radius. Line strength are consistent with Sanc hez-Blazquezl 
(2004). 

NGC4150: Noticeable obscur ation from dus t is visible in the 
central parts of this galaxy l Ouillen et al. 2000). A kinematically 
decoupled component resides in the central few arcseconds of this 
galaxy (see Paper III), coincident with a strong drop in the Mg b 
absorption strength and a corresponding peak in H/3 absorption 
strength, which itself is globally enhanced. The iron line maps are 
relatively featureless by comparison, with no strong central fea- 
tures. The core is therefore domina ted by a young stellar popula- 
tion, as found by other authors (e.g., Fernan das et all2004l) . 

NGC4262: This strongly barred object shows inconspicuous 
metal line gradients. The H/3 absorption map shows a mild increase 
with radius. 

NGC 4270: The central peaked metal absorption line strengths in 
this galaxy show evidence of an increase beyond 10" along the ma- 
jor axis, giving a small dip around 4-5" from the center. The H/3 
absorption map shows no significant features. 

NG C 4278: This LINE R <Ho et alj 1 19971) and radio source 
(e.g.. ICa petti et a f1l2000t) shows some central dust features (e.g., 
IPeng et alj|2002l) . There is some evidence that the Mgfe contours 
are more flattened than the isophotes. There is a well-defined de- 
pression in the H/3 line strength within the central 5" (see also 
iDavies et alJI 19931) which we partly ascribe to imperfect removal 
of strong H/3 emission. 

NGC 4374 (M84): This well known g iant elliptical shows notice- 
able dust absorption l Ouillen et al. 2000) and has a BL Lac nucleus 
iBower et alj|200ol) . The metal lines have a similar morphology to 
those in NGC 4278, following the isophotes in general. The H/3 ab- 
sorption map is rather flat, showing a slight increase towa rds the 
edge of the field, a lso p resent in the data o f [Day ies et alj ll993l) . 
ICarollo et afl<1993l) and lSanchez-BlazquezI fiwfi) . 

NGC 4382 (M85): This well known object shows a distinct de- 
pressio n in Mg b within the central 10", consistent with lFisher et alj 
1 1996), which appears slightly more flattened than the isophotes. 
This is coincident with an enhancement in the H/3 absorption map 
(see also lKuntschner et all200l|) . as well as corresponding closely 
to a region of peculiar kinematic behavior (decoupled rotation, 
low velocity dispersion, and complex higher-order terms - see Pa- 
per III). There is also evidence for a flattened feature in the iron ab- 
sorption maps within this region, most evident in the Fe5015 map. 

NGC 4387: This boxy galaxy iPeletier et alll990l) shows regular 
metal line gradients consistent with the isophotes. The H/3 map ap- 
pears relatively flat in the inner regions while we observe a drop of 
H/3 strength towards larger radii. 

NGC 4458: This galaxy with a kinematically decoupled compo- 
nent (hereafter, KDC; see Paper III) shows line strength gradients 
consistent with the isophotes . The region of the K DC is inconspic- 
uous in the H/3 map (see also Morelli et al. 2004). 

NGC 4459: This galaxy shows a central region of stronger H/3 
absorption with a ring like structure around it. This is consistent 
with the central dust ring found in this galaxy; see Paper V and 
ISarzi et ail 1200 ll) . The metal lines show relatively regular gradi- 
ents. 

NGC 4473: This galaxy exhibits a complex morphology in the 
velocity dispersion map, with a region of high dispersion along 
the major-axis which widens at larger radii, which probably corre- 
sponds to two counter-rotating stellar disks. The metal line strength 



maps are more flattened than the isophotes, indicating a higher 
metallicity in the disk(s) while the H/3 map is rather flat. 
NGC 4477: This galaxy shows a prominent misalignment of the 
kinematic and photometric major-axis. The metal line strength 
maps show an inner region with enhanced metallicity and show lit- 
tle structure elsewhere. There is weak evidence at larger radii along 
the major axis for increased metal line strength. The H/3 map ap- 
pears flat. 

NGC 4486: The line strength in the central 3" and some part of 
the jet region of this famous Virgo galaxy are clearly affected by 
an imperfect removal of emission lines and non-thermal emission 
iDavies et alll993l) . Overall, the line strength maps appear flat over 
large parts of the FoV. The H/3 map exhibits notably low values. 

NGC 4526: This SO object has a prominent dust disc which is not 
only visible in the reconstructed SAURON image but also clearly 
seen in the Mg b and H/3 maps. The strong H/3 line strengths in the 
region of the dust lane indicates the presence of dust enshrouded 
young stars. 

NGC 4546: This galaxy exhibits ionized gas counter-rotating 
with respect to the stars (see Paper V). The Mg b contours are more 
elongated in the direction of the rotation than the isophotes. This 
effect is less prominent in the iron maps. The H/3 map appears rel- 
atively flat with a hint of weaker H/3 absorption along the major 
axis. 

NGC 4550: This SO galaxy has two counter-rotating stellar disks 
of similar mass iRix et al. 1992). There is evidence for the Mgfo 
contours being flatter than the isophotes. The H/3 map shows a 
mildly elevated level in the central region along the major axis. 

NGC 4552: This large elliptical shows regular line strength gradi- 
ents, c onsistent with lSanchez-BlazouezI 120041) and Rampazz o et alj 
1 2005). The H/3 map appears flat and shows notably low values. 

NGC 4564: This highly elongated galaxy shows isoindex con- 
tours of the Mg b and the Fe indices which are more flattened than 
the isophotes, indicating that the bulge and disk of this galaxy have 
different stellar populations. 

NGC 4570: This edge-on SO galaxy features multiple kine- 
matic com ponents s uch as an outer and a nuclear disk 
Ivan den Bosch & Emsellemlll998l) . The Mg b map and to a lesser 
extent the Fe maps show a central peak surrounded by a region de- 
pressed in line strength. At larger radii the onset of the outer disk, 
showing strong metal lines, can be seen. The H/3 map is relatively 
flat. 

NGC 4621: This is another object with a strong indication from 
the stellar kinematics of a disc component, which is also seen in the 
metal line strength maps showing more flattened isoindex contours 
than the isop hotes. The inner 60 pc harbor a counter-rotating core 
iWernli et al. 2002). The H/3 map appears relatively flat. 

NGC 4660: Similar to NGC 4621, there are strong indications 
for a disk component in this galaxy, accompanied by pointy disky 
isophotes. This galaxy is among the best examples of Mg b isoin- 
dex contours to be flatter than the isophotes for galaxies with strong 
rotation. 

NGC 5198: This galaxy, with a central KDC, rotating nearly per- 
pendicularly to the outer body, shows regular metal line gradients 
and a mild positive H/3 gradient. 

NGC 5308: This is a typical case of a close to edge-on disc galaxy 
showing a rapidly rotating component (Paper III) and a double 
sign r eversal in /13, indicating a central bar iBureau & Athanassoulal 
2005). The metal line strength gradients are regular and appear to 
be consistent with the isophotes. 

NGC 5813: A galaxy with a well-known KDC jEfstathiou et alj 
ll98lll98llBender et aiTl9 94. Paper III). The line strength maps 
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of this galaxy have already been presented in Paper II. Here it 
is shown tha t we re produce the large Mgfe gradient presented by 
iGorgas et alJ ll990l) . The KDC is not seen in the absorption line 
maps. 

NGC5831: The ki nematics of this galaxy reveal the p resence of 
a well-known KDC iDavies etal]|l983t IPeletier et al1ll99(l Pa- 
per III). The KDC is not really seen in the line strength maps, apart 
from the Mg b map which appears to be more elongated in the direc- 
tion of the KDC than the isophotes. The H/3 map appears relatively 
flat. 

NGC 5838: This boxy galaxy is a strong rotator and lPeletier et alJ 
1 1990) show that the bulge of this galaxy contains old stellar pop- 
ulations with small color gradients, except for the inner 2". This is 
consistent with our line strength maps which show regular gradients 
and a mild peak of the H/3 map in the very center. There is evidence 
for the Mg b contours being flatter than the isophotes. Our emission 
line maps (see Paper V) also show significant detection in the center 
with evidence for a rotating structure. 

NGC 5845: This galaxy is a compact elliptical, close to 
NGC 5846. S mall gradients in the line strength maps are seen, 
consistent with Sanchez-Blazquez (2004). The central disk, clearly 
seen in the stellar kinematics (Paper III), is not apparent in the line 
strength maps. 

NGC 5846: The SAURON field of view of this bright giant ellipti- 
cal includes a foreground star and a companion (NGC 5846A) north 
and south of its nucleus, respectively. The line strength maps are 
generally regular , with small gradient s. T he line strengths are con - 
si stent with Sanc hez-Blazquezl J2004 ) and Rampazzo et alji2005l) . 

NGC 5982: F or this galaxy we co nfirm the presence of a KDC 
first detected by Wa gner et alJ<1988l) . The KDC is not clearly seen 
in the line strengths maps, however, there is evidence for the Mg b 
contours being flatter than the isophotes. The H/3 map appears rel- 
atively flat. 

NGC 7332: Another strongly bo xy galaxy with a KDC (c ounter- 
rotating) in the central 3" (see Falco n-Barroso et"al] Ebo4). There 
is evidence for the Mgfe contours being flatter than the isophotes. 
The galaxy shows high H/3 values, indicating an age of about 
3 Gyr (se elFalcon-Barroso et alj20 04f). The data are co nsistent with 
ISanchez-BlazaueJ 120 04) and Rampa zzo et alJ J2005). 

NGC 7457: This small SO galaxy with a KDC shows surprisingly 
little struc ture in its line strength maps, consistent with the color 
profiles of IPeletier etal] jl999). The high H/3 values indicate a rel- 
atively young age. 
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